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Abstract 

The curvature stress of DNA packed inside a phage is balanced against its 
electrostatic self-interaction. The DNA density is supposed nonuniform and as a 
result the Donnan effect is also inhomogeneous. The coarse-grained DNA density 
is a nonlinear function of the DNA radius of curvature at a given position inside 
the bacteriophage. It turns out that a region (or regions) exists totally free from 
DNA. The size of such holes is computed. 
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The packing of DNA in cells and viruses are problems in the nanothermodynamics 

of polyelectrolytes i.e. the systems are small and hence the free energies need not be 
extensive. A case in point is the compaction of DNA inside bacteriophages where the 
extensive self-energy of the DNA spool is balanced against the non-extensive energy of 
curvature arising mainly from a region of tightly wound DNA surrounding a hole [2]. 
The degree of non- extensivity depends on the problem at hand. For fairly large DNA 
globules that are compacted in solution, the shape is determined by non-extensive terms 
whereas the volume is not j^l E] . 

In previous work [21 El E] we assumed the DNA density in a coarse-grained continuum 
approximation to be homogeneous, but we relax this supposition here. This introduces 
a peculiar effect : the Donnan effect is now also nonuniform and it happens to compete 
with the curvature stress which increases sharply in tightly wound regions. Although the 
uniform approximation of [2| agrees well with the DNA spacings within bacteriophage T7 
determined experimentally 0, the forces encountered in single molecule experiments jH] 
of DNA being forced into Phi 29 phage by the connector motor seem to overestimate the 
predicted forces [21 U\ ■ Thus a reinvestigation of at least one assumption is warranted. 

Close-packed DNA in vitro has often been found to be a hexagonal liquid crystal [2111]. 
Sometimes, for long polydisperse DNA, the hexagonal phase appears to be precluded in 
favor of a possible "hexatic" phase JUj but a full investigation of the algebraic decay of the 
ovientational order has yet to be carried out. Here, for the sake of definiteness, we assume 
the DNA lattice inside a phage is hexagonal (for general, biologically oriented reviews on 
DNA packaging in bacteriophages, fn\ IT^). However, such a lattice is probably rarely 
ever ideal but is bound to contain DNA crossovers leading to knotted entanglements upon 
release of the DNA from the phage ITlj . 

Theories of slightly deformed hexagonal phases have been set up previously I16j . 
They apply to macroscopic or extensive systems. Here, our compacted DNA is of meso- 
scopic dimensions and its properties are dominated by nonextensive curvature stresses 
so we simply disregard fluctuations in the director n{r) and coarse-grained DNA density 
p(r) at position r. The linear double-stranded DNA has only two ends so to a good 
approximation a well-known continuity law applies [T7| I18j. 

V-pn = (1) 
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We next assume absence of splay ( V ■ n = 0) since the splay elastic constant is expected to 
be very large. In general, at constant density, we would deal with the theory of so-called 
developable domains in describing the deformed hexagonal liquid crystal ll^. We here 
consider deformations of nonuniform densities in such a way that the director remains 
perpendicular to the density gradient. Under these conditions eq.(IH) is obeyed. 

Now the simplest total free energy of the DNA coil as a functional of the density p(r) 
seems to be 

Ftot Fint 

W)] + U, W)] (2) 

where Fi^t is the free energy of the DNA interacting with itself in a hypothetically straight- 
ened hexagonal lattice (though it is inhomogeneous across the plane perpendicular to its 
central axis) and Uh is the bending energy given by 

U, = ]-PkBT ("^ ds R:\s) (3) 

Z JO 

Here, P is the DNA persistence length, ks is Boltzmann's constant, T is the temperature 
and Rc is the radius of curvature of the DNA wormlike coil at position s on the DNA 
contour of total length L. As already stated, undulations of the DNA will be disregarded 
in the tightly packed state we focus on. If we define the surface area S{t) of the unit 
cell of the deformed hexagonal lattice at position r, we may introduce a continuum 
approximation ^(r)^^ = dr and p(r) = l/AS{v). The length of a DNA nonomer is A so 
that p(r) is a monomer density from now on. 

U, = \pAksT ( dv-^ (4) 

The number of monomers inside the phage of volume V is kept fixed. 

N = ^ = f drp(r) (5) 
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Supposing gradients of p(r) may be disregarded in Fint, we now minimize the total free 
energy with respect to p(r) while accounting for the constraint eq.(j^. 

, , PAkeT 

^^^^ ^ 2R'^{t) ^ ^ ^ ' 

We have introduced the inhomogeneous bulk chemical potential of the DNA monomers 
(p = 6Fint/6p{r)). It is readily seen that the curvature term becomes appreciable for 
tightly wound regions appropriate for DNA packed in phages. 



Next, we need an expression for the DNA chemical potential. The phage is in thermo- 
dynamic equilibrium with a large reservoir containing a univalent salt of concentration 

Cs- 

The phage is presumed to be free from multivalent ions. We first focus on a ho- 
mogeneous hexagonal lattice of DNA and its Donnan equilibrium with such a reservoir. 
One convenient approximation for the osmotic pressure is due to Oosawa pUl I21j . For 
our purposes, we set the osmotic coefficient equal to a constant so we write the osmotic 
pressure tTos of the uniform bulk hexagonal phase with respect to the buffer as 

/ ApokBT\ 



h{w) (7) 



w = 2.?>^HlQcs = ASoQcs 

h{x) = {I + x^)^ - X (8) 

The uniform monomer density of the lattice is po and the Bjerrum length Q = q^/eksT 
in terms of the elementary charge q and the permittivity of the solvent e. For a uniform 
hexagonal lattice, the area of the unit cell is 5*0 = 3^Hq/2 where Hq is the lattice spacing 
i.e. the distance between centeraxes of adjacent DNA helices (po = iM'S'o). Eq.((7j) is, of 
course, an approximation though, on the whole, it does agree with recent careful mea- 
surements of the osmotic pressure of DNA suspensions by Raspaud et al [221 (to a good 
approximation, the pressure conforms to a universal form tTos = pokBToj^Cg/ po) indepen- 
dent of any model; the osmotic coefficient is also quite constant). At very high concen- 
trations, the Oosawa approximation does seem to underestimate experimental pressures 
This has been attributed to a failure of the Oosawa model at close separations [21] 
but a calculation not invoking cylindrical symmetry still needs to be presented. Under 
conditions of tight packing, it is reasonable to neglect undulatory entropy. The bulk free 
energy is now the work of compressing the lattice from a state of infinite separation to 
the state at hand all the while keeping the chemical potential of the reservoir electrolyte 
fixed. 

Pint [po] = Lfint + const.N = -3H r dH'H'TTosiH') + const.N (9) 

J oo 

Unt = (^]nw) (10) 
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Hence, the chemical potential /i = dFint [po] / dN of a DNA monomer at a fixed volume 
and buffer concentration may be expressed by 



AksT ^ 
const. H -— m 

4g 



:i+w2)5 - 1 



(12) 



Note that the limiting forms of the logarithmic term are —2lnw as -C 1 and 2/w as 
w ^ 1] still , in practice, the variable w is often of order unity. 

We next apply eq.lfT^ to the deformed DNA lattice inside the phage whose density 
p(r) = 4:Qcs/Aw{r) at equilibrium is given by eq.®- 



;i + «;2(r))5 + l 



(l + u;2(r))i - 1 
where A; is a constant. This is rewritten as 

'2Qc 
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(13) 



(14) 



A ) Zt(r) 

implying that the dimensionless function Z{j) must be larger than unity otherwise a 
solution to eq.(^ni) does not exist. Without knowing the precise structure of the DNA 
globule inside the phage, we already conclude from eq.lfT^ that tightly wound regions of 



curvature radius Rc — \f2PQ ~ 8 nm must be depleted of DNA (persistence length P = 
50 nm; Bjerrum length Q = 0.71 nm). The tight bending of the DNA causes an increase 
in the curvature stress which has to be relieved by lowering the electrostatic stress (see 
fig. 1). 

Let us now consider the profile about a cylindrical hole within a phage filled with DNA 
of genomic size. Then, the DNA almost fills the compartment of volume V completely 
so we have from eq.© 

T , , , 2QcsV{k - 1) 

dr p(r) ^ 4 (15) 
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By way of example, let us reconsider the data on phage T7 of Cerritelli at al j3] which we 
previously analyzed in terms of a model of uniform density [2]. The NaCl concentration 



of the buffer was 0.1 M. The DNA genome of length L is enclosed in a volume V of about 
77600 nm^ which almost equals the volume V^na of the DNA lattice itself. Eq. (fTH|) then 
yields X = 2.06 and k = 6.08. Accordingly, the radius E of the (presumably single) hole 
free from DNA is computed from eq. lfT^ : 

yielding ~ 6 nm. This should be compared with a hole radius of about 4 nm or 
possibly larger, based on the experimental spacing {H = 2.54 nm) and estimates of the 
inner volume of bacteriophage T7 [21 E] • An alternative route to a theoretical estimate 
for E is via the (average) experimental spacing H: the volume Vdma of the DNA lattice 
is slightly less than V and so ~ LS. Hence, we have X = 2.10 from eq. ffTBj) and again 
^ 6 nm. Similar conclusions apply to recent structural work on bacteriophage T4 [25]. 
If the phage is modeled as an ideal sphere of radius i? = 41.5 nm, the radius of the hole 
would be about 10 nm (setting L = 57800 nm; Cg ~ 0.3 M) from eq. (fT7j) . Olson et al 
studied the DNA density which decreases inward starting from the outer capsid shell and 
there does appear to be a hole of radius ~ 10 nm. However, at present, it is difficult to 
make this comparison more quantitative. Finally, we can justify the neglect of density 
gradients and Coulomb integral terms in eq.©; they are only of relative magnitude 
H'^/Rl which is much smaller than unity. 

We conclude from a simple model of close-packed DNA enclosed in phages, in which 
we relax the supposition of uniformity, that a hole (or holes) must occur free from DNA, 
ionic strengths typical of laboratory or physiological conditions. The curvature 
stress in a tightly wound region of DNA becomes so high that the nonuniform Donnan 
effect can only compete up to a point; the DNA is thus forced out of such a region beyond 
a certain curvature given by eq. lfTTj) (See fig. 2). It would be interesting to see if the 
DNA density may be approximated by eq. p3|) - which is independent of the shape of 
the capsid (!) - but one has to be wary about the limitations of the present continuum 
analysis. At low densities, toward a hole, the DNA mesophase should become more 
disordered. Kasspidou and Van der Maarel have established that a hexagonal phase of 
short DNA rods melts into a cholesteric phase jH]. This happens, for instance, at an 
interaxial spacing oi H = 3.7 nm in a 0.10 M NaCl solution. It is then predicted from 
eq. lfT^ that such a transition would occur at a radius of curvature i?c = 9 nm in the 
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T7 phage [Sj, at least if we are allowed to extrapolate the results for short rods P to 
long DNA. At such nanoscales the concept of phases becomes very tentative of course. 
Anyway, the more disordered cholesteric "phase" surrounding the hole would still have 
an average spacing H very similar to that of the coexisting hexagonal phase in the 
bulk of the phage. Another possibility for disorder is an escape of a section of the DNA 
"into the third dimension" i.e. if a DNA strand aligns along the cylindrical axis of the 
hole. Such a state would seem to be unstable however. 

The buffers used in jHl I2S1 also contain magnesium ions though not in excess. When 
genomic DNA is then packaged in bacteriophage T7, two Raman bands of the DNA are 
significantly perturbed with respect to DNA in solution |2n]- Overman et al attribute 
this to enhanced magnesium concentrations within the phage j^H]- The effect of divalent 
ions has not been allowed for here although we note that magnesium is absent in the 
buffers used by Cerritelli et al 0. 

Small angle X-ray scattering could be useful in trying to verify the profile predicted 
by eq. ()14|) . However, it is well known that the scattering by phages is subtle with regard 
to interpretation j27j. The wavevector dependence of the scattering by lyotropic liquid 
crystals is peculiar in its dependence on the osmotic compressibility and the elastic moduli 
|28j . It might be useful to extend theories pertaining to scattering from the bulk ^3 
IT^ ITH] to that from the DNA mesophase confined in phages. We remark that there is a 
length scale in the theory for a fluctuating splayless liquid crystal which turns out to be 
identical to our decay length ^JIPQ in eq. lfT^ (in [131 CSl CHI , there is also a competition 
of bending and osmotic stresses). Therefore, the DNA liquid crystal near the outer rim 
of the phage behaves much like a useful polymer liquid crystal but this is no longer true 
near the hole. 
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Figure 1: Schematic representation of DNA packed in a spherically shaped virus (plane 
perpendicular to the spool axis) . The packing is locally hexagonal and becomes less close 
in the inner region of tight winding. 



10 



p 




E r 

Figure 2: Coarse-grained DNA density in a plane through the central axis (MO) of a 
hole within a presumably spoollike region of DNA inside the phage. The size of the 
nonuniform region beyond E is approximately equal to ^2PQ. 
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